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Abstract 

Parallel flow in a Hele-Shaw cell occurs when two immiscible liquids flow 
with relative velocity parallel to the interface between them. The interface 
is unstable due to a Kelvin-Helmholtz type of instability in which fluid flow 
couples with inertial effects to cause an initial small perturbation to grow. 
Large amplitude disturbances form stable solitons. We consider the effects of 
applied magnetic fields when one of the two fluids is a ferrofluid. The disper- 
sion relation governing mode growth is modified so that the magnetic field 
can destabilize the interface even in the absence of inertial effects. However, 
the magnetic field does not affect the speed of wave propogation for a given 
wavenumber. We note that the magnetic field creates an effective interaction 
between the solitons. 
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The Saffman- Taylor problem |^ considers two immiscible viscous fluids moving in a 
narrow space between two parallel plates (the so-called Hele-Shaw cell). When a low viscosity 
fluid invades a region flUed with high viscosity fluid, the initially flat fluid-fluid interface is 
unstable and evolves through a mechanism known as viscous flngering [||. We call the 
displacement of one fluid by another frontal flow. In contrast, parallel flow occurs when the 
fluids flow parallel to the interface separating them. One important example of parallel flow 
occurs after the passage of a fully developed Saffman- Taylor flnger. 

Recent experimental and theoretical studies examined the dynamics of fluid in- 

terfaces under parallel flow in Hele-Shaw cells. Zeybek and Yortsos studied, both 
theoretically and experimentally, parallel flow in a horizontal Hele-Shaw cell in the large 
capillary number limit. For flnite capillary number and wavelength, linear stability analy- 
sis indicates that small perturbations decay, but the rate of decay vanished in the limit of 
large capillary numbers and large wavelength. Furthermore, a weakly nonlinear analysis of 
the problem found Korteweg-de Vries (KdV) dynamics leading to stable flnite amplitude 
soliton solutions. Solitons were indeed observed experimentally. Gondret and Rabaud 
incorporated inertial terms into the equation of motion in a Hele-Shaw cell and found a 
Kelvin-Helmholtz instability for inviscid fluids. For viscous fluids they derived a Kelvin- 
Helmholtz-Darcy equation and found the threshold for instability was governed by inertial 
effects, while the wave velocity was governed by the Darcy's law flow of viscous fluids. Their 
experimental results supported their theoretical analysis. 

As was the case for frontal flow of nonmagnetic fluids in Hele-Shaw cells, many research 
groups have studied the frontal interface behavior when one of the fluids is a ferrofluid 0, 
and an external magnetic fleld is applied 0-lO|. Ferrofluids, which are colloidal suspensions 
of microscopic permanent magnets, respond paramagnetically to applied fields. As a result 
of the ferrofluid interaction with the external fleld, the usual frontal displacement viscous 
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fingering instability is supplemented by a magnetic fluid instability , resulting in a variety 
of new interfacial behaviors. Depending on the applied field direction, one observes highly 
branched, labyrinthine structures patterns showing an ordered line of peaks ||^, or even 
the supression of the usual viscous fingering instability |jTO|. Rosensweig discusses the 
Kelvin-Helmholtz instability for unconfined ferrofiuids. 

In this paper we perform the linear stability analysis for parallel fiow in which one fiuid is 
a ferrofiuid and a magnetic field is applied. We consider three separate field configurations: 
(a) tangential, for in-plane fields tangent to the unperturbed interface; (b) normal, for in- 
plane applied fields normal to the unperturbed interface; (c) perpendicular, when the field is 
perpendicular to the plane defined by the Hele-Shaw cell plates. We show the magnetic field 
provides additional mechanisms for destabilizing the interface, and we analyze qualitatively 
the interactions between solitons caused by the magnetic field. We neglect inertial terms 
because they are not needed to understand the interfacial instability. 

Let us briefiy describe the physical system of interest. Consider two semi-infinite immis- 
cible viscous fiuids, fiowing with velocities Ui and U2, along the x direction, in a Hele-Shaw 
cell of thickness b (see figure 1). We assume that b is smaller than any other length scale 
in the problem, and therefore the system is considered to be effectively two-dimensional. 
Denote the densities and viscosities of the lower and upper fiuids, respectively as pi, rji 
and P2, r]2- To achieve steady-state parallel fiow the velocities and viscosities must obey 
the condition rjiUi = rj2U2- According to Gondret and Rabaud 0, we may neglect inertial 
terms relative to viscous terms provided kb « 12/Re, where is a typical wavevector and 
Re = pUb/rj is a characteristic Reynold's number. 

Between the two fiuids there exists a surface tension a. We assume that the lower fiuid is 
the ferrofiuid (magnetization M), while the upper fiuid is nonmagnetic. In order to include 
the acceleration of gravity g, we tilt the cell so that the y axis lies at angle (3 from the 
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vertical direction. To include magnetic forces, we apply a uniform magnetic field which 
may point along the x, y or z axis. During the flow, the fluid-fluid interface has a perturbed 
shape described as y = C{x,t) (sohd curve in figure 1). 

Hydrodynamics of ferrofluids departs from the usual Navier-Stokes equations through the 
inclusion of a term representing magnetic force. Let M represent the local magnetization 
of the ferrofluid, and note that the force on M depends on the gradient of local magnetic 
field H. The local field differs from the applied field Hq by the demagnetizing field of the 
polarized ferrofluid. We will assume M takes a constant value parallel to the applied field. 
This amounts to neglecting the demagnetizing field relative to the applied field and can 
be justified for low magnetic susceptibility of the ferrofluid, or for large applied fields that 
saturate the ferrofluid magnetization. It can also be justified for very thin ferrofluid films 
when the field is parallel to the plane of the cell. 

For the quasi two-dimensional geometry of a Hele-Shaw cell, the three dimensional flow 
may be replaced with an equivalent two-dimensional flow v{x, y) by averaging over the z 
direction perpendicular to the plane of the Hele-Shaw cell. Imposing no-slip boundary 
conditions and a parabolic velocity profile one derives Darcy's law for ferrofluids in a Hele- 
Shaw cell OO, 



= -j^ I Vp - - J_^^^ (M ■ V)Hdz - p{g ■y)y\, (1) 

where p is the hydrodynamic pressure. Equation (0) describes nonmagnetic fluids by simply 
dropping the terms involving magnetization. 

When the velocity field v is irrotational, it is convenient to rewrite equation (|l|) in terms 
of velocity potentials. We write v = — V0, where defines the velocity potential. Likewise 
we introduce the scalar magnetic potential 

if = / T^^d\' 2 
Js r — r 



where H = — Vy^. Here the unprimed coordinates f denote arbitrary points in space. The 
primed coordinates f are integration variables within the magnetic domain S, and ci^r' 
denotes the infinitesimal area element. The vector n' represents the unit normal to the 
magnetic domain in consideration. 

To study the interface dynamics, we evaluate equation (|I|) for each of the fluids on the 
interface, subtract the resulting equations from each other, and divide by the sum of the 
two fluids' viscosities to get the equation of motion 



2 ) \ 2 ) 12(771 + 772) 

1 r+Va ^ ^ ] 
crK + - {M ■Vip)dz + {p2 - pi)g cos P y\. 

J-b/2 J 
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To obtain (j^) we have used the pressure boundary condition P2 ~ Pi = crn at the interface, 
where k = {d'^(/dx'^)[l + (SC/Sx)^]"^/^ denotes the interfacial curvature in the plane of 
the Hele-Shaw cell. The dimensionless parameter A = (772 — f]i)/{i]2 + f]i) is the viscosity 
contrast. 

We perturb the interface with a single Fourier mode 

C{x,t) = Cfc exp {i (ujt-kx)). (4) 

The velocity potential for fluid z, (^j, must contain the uniform unperturbed flow Ui and a 
perturbed part that reflects the space and time dependence of C, obeys Laplace's equation 
V^0j = and vanishes as ?/ ^ ±00. The velocity potentials obeying these requirements are 



exp(±|fc|?/) exp(2(a;t — kx)) — Ui x. (5) 



To conclude our derivation and close equation (^) we need additional relations expressing 
the velocity potentials in terms of the perturbation amplitudes. To flnd these, we considered 
the kinematic boundary condition, which states that the normal components of each fluid's 



velocity at the interface equals the normal velocity of the interface itself Inserting 
expression (|) for C{x,t) and (|D for 0i into the kinematic boundary condition, we solved 
for (f)ik{t) consistently to first order in ( to find 

iuCk k 

= — 1^ + * 1^ ^iCfc) (6) 

and a similar expression for (j)2k- 

Substitute expression (||) for (pik and the related expression for 02fc into equation of 
motion and again keep only linear terms in the perturbation amplitude. This procedure 
eliminates the velocity potentials from equation (|^), and we obtain the dispersion relation 
for growth of the perturbation ({x, t) 
f ViUi + r]2U2\ 



UJ 



NBl,{k) - {kbY - {koby (7) 



V V1+V2 J 12(r7i + 
where A''^ = 2M%/a is the magnetic Bond number and fco = y [(pi — P2)g cos f3]/a. 



The real part of u; is /c times the phase velocity, and is the viscosity-weighted average 
of the two fluid velocities. Note that the magnetic field does not alter the phase velocity 
of the waves. The imaginary part of u, which governs the exponential growth or decay 
of the wave amplitude, does include effects of the magnetic field. Exponential (unstable) 
growth occurs when the imaginary part of u is negative. We point out that when there is no 
applied magnetic field (A''^ = 0) our equation (|^ agrees with the dispersion relation derived 
by Gondret and Rabaud for the case in which the cell is vertical [jS = 0) and by Zeybek 
and Yortsos [^,0 for the case in which the cell is horizontal (/? = vr/2). 

Terms containing Ij{k) originate from the Fourier transforms of 

the magnetic contribution to equation (H). The subscript j = x,y, z indicates the tangential, 
normal and perpendicular magnetic field configurations, respectively. For M in the x or y 
direction we can expand equation (H) to first order in ( to obtain 
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and 



dx'{x — x') 



dC{x') 
dx' 



F{x - x') 



(9) 



Iy{x) 



dx'lCix) - Cix')]F{x - x') 



(10) 



where 



Fix) = 



I i'+b/2 j-+b/2 



dzdz' 



b J-b/2 J-b/2 [x^ + (-2 — -Z')^]^/^ bx' 



\x\ 



(11) 



In contrast, for Iz{x) the z integration inverts the derivative of (f with respect to z in 
equation (H) so that after integrating over y' and expanding to first order in powers of (, 
this term simphfies to 



Fix) 



oo 2 

dx' - 

-oo b 



(12) 



\x — x'Y y (x — x'y + b'^ 
We obtain the specific forms for the Jj(/c)'s corresponding to each particular field config- 
uration by taking the Fourier transform of equations (PD]), and ([T^). After some simple 
algebra we find the following expressions for the magnetic terms Ij{k) 



(^) [^{kby + r2 - r] rfr, 



(13) 



(14) 



and 



h{k) = 4 / sin'r 
Jo 



dr. 



J ^{kb/2f + r2 

In the limits of small and large wavevector these Fourier transforms reduce to 



(15) 



h{k) 



-[(3/2-C + ln2) -lnA;6](A;6)2 kb«l 



(16) 



—Ttkb 



kb » 1 



[(2-C + ln2) -lnA;6](A;6)V2 kb « 1 
Iyik)^{ (17) 

nkb kb » 1 



\{l-C + \n2)-\nkb]{kb)y2 kb « 1 
h{k) ~ <; (18) 

\n{kb/2) kb»l, 



where C ~ 0.57721 denotes Euler's constant |T3|. Our results (pTSD , (|l^ and ([T5|) agree 
with similar kind of calculations related to frontal displacements in Hele-Shaw cell with 
ferrofluids 0-|lO|. 



The dispersion relation (0) is given for the case of systems with infinite extent along the 
y-axis. For finite extent L the algebraic dependence on wavevector k is modified by a first 
order rational function of sinh kL as shown by Zeybeck and Yortsos . When kL is large 
this finite size correction dies off exponentially quickly. The magnetic integrals Ij{k) likewise 
possess exponentially small finite size corrections. 

Consider the stability of the fluid-fluid interface for the different field configurations. The 
initially fiat interface is unstable to perturbations with wavenumber k when NBlj{k) — {kbY — 
{koby is positive. If the heavier fluid is below the lighter fluid, (pi > P2), then both gravity 
and surface tension stabilize the system and ko is real. Therefore, in the absence of applied 
magnetic field {Nb = 0), the temporal growth rate of any perturbation is negative and waves 
are damped. On the other hand, if the external magnetic field is nonzero, the stability of 
the interface will depend on the field's direction. Figure 2 illustrates how the magnetic 



terms (0), (0) and ( [Tsl) vary with reduced wave number kb. Inspecting figure 2 and 
the imaginary part of the dispersion relation (0) we note that a tangent field configuration 
{Ix{k) < 0), makes the growth rate even more negative than when the field is absent. So a 
tangent external field has a stabilizing nature, reinforcing the effects of gravity and surface 
tension. In contrast, since Iy{k) and Iz{k) are both positive quantities, if a sufficiently 
strong magnetic field is applied normal to the fiuid-fiuid interface, or perpendicular to the 
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cell plates, the growth rate may become positive, leading to a possible destabilization of 
the interface. We conclude that the magnetic field can destabilize the interface even in the 
absence of inertial effects. 

In addition to the interface stability issue discussed above, it is interesting to ask how the 
magnetic field acts on the motion of interfacial waves once they appear. In the following, we 
discuss the action of the applied magnetic field on the solitons that appear in parallel flow 
in Hele-Shaw cells. To treat the problem rigorously would require reproducing the analysis 
of Zeybek and Yortsos that derived Airy and KdV equations from a weakly nonlinear 
analysis of the interfacial perturbations. Here we simply point out that the solitons may be 
considered as localized perturbations on the fiat interface. When magnetic fields are present 
the solitons acquire net dipole moments equal to the magnetization of the fiuid multiplied 
by the integrated area of the soliton. 

Take the generic form of a KdV soliton, 

t) = -^sech (^^{x - ct) j , (19) 

written here in terms of the scaled time, position and height variables discussed in [^], where 
c is the speed of propogation. We define the scaled dipole moment of the soliton of speed c 
as 

/oo _, 
Mu{x, t)dx = - VcvrM. (20) 
-oo 

In doing so, we neglect the magnetic field dependence of the shape of the soliton. We may 
consider the magnetic moment (PDj ) as the leading, linear term in a perturbative series in 
powers of applied field, and expect a cubic correction due to the field-dependent soliton 
shape. As noted in [Q, the actual profile in unsealed coordinates may be either positive 
or negative, and the dipole moment given here must be divided by the position and height 
rescaling factors to yield the true moment. True dipole moments m point parallel to the 

9 



magnetization M when the sohton consists of excess magnetic fluid, and points opposite to 
M when the sohton consists of missing magnetic fluid. 

Dipole interactions are long-ranged, falling off as for moments separated by a dis- 
tance X. This contrasts with the fluid-dynamic interaction of solitons which decays expo- 
nentially with separation. An interesting additional feature of the dipole-dipole interaction 
is its variation with the relative orientation of dipole moments and the vector joining them. 
In the case of solitons with parallel moments rhi and m2 displaced from each other along 
the X axis, the interactions will be attracting, with strength 2mim2, when the magnetiza- 
tions lie along the x axis (tangetial) but will be repelling, with strength mim2 when the 
magnetizations lie along the y (normal) or z (perpendicular) axes. 

In conclusion, we have performed the linear stability analysis for parallel flow in a Hele- 
Shaw cell when one of the fluids is a ferrofluid. We show that the magnetic held may 
provide a new mechanism for destabilizing the interface in the absence of inertial effects, 
and we determine the magnetic correction to the dispersion relations for three distinct fleld 
orientations. Finally, we suggest parallel flow of ferrofluids as a novel system in which to 
investigate soliton interactions. 
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FIGURE CAPTIONS 



FIG. 1: Schematic configuration of the parallel flow geometry. 

FIG. 2: Variation of Ij{k) as a function of kb for (a) tangential, (b) normal, and (c) perpen- 
dicular magnetic field configurations. 
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